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We determine, using a self consistent method, the charge and spin currents in ballistic Josephson junctions 
consisting of several ferromagnetic (F) layers sandwiched between superconducting (S ) electrodes (S FS -type 
junctions). When there are two F layers, we also consider the experimentally relevant configuration where a 
normal (N) nonmagnetic spacer layer separates them. We study the current-phase relationships as functions 
of geometrical parameters that are accessible experimentally including particularly the angles that characterize 
the relative orientation of the magnetization in the F layers. Our self-consistent method ensures that the proper 
charge conservation laws are satisfied, and that important proximity effects are fully and properly accounted 
for. We find that as we vary the phase difference A(p between the two outer S electrodes, multiple harmonics 
in the current phase relations emerge, the extent of which depends on the interface scattering strength and on 
the relative F layer widths and magnetization orientations. By manipulating the relative F layer magnetization 
orientations, we find that the charge supercurrent can reverse directions or vanish altogether. These findings are 
discussed in the context of the generation and long-range nature of triplet pair correlation within these structures. 
We also investigate the spin currents and associated spin transfer torques throughout the entire junction regions. 
For noncollinear relative magnetizations, the non-conserved spin currents in a given F region gives rise to net 
torques that can switch directions at particular magnetic configurations or Aip values. The details of the spin 
current behavior are shown to depend strongly on the degree of magnetic inhomogeneity in the system, including 
the number of F layers and the relative widths of the F and N layers. 


I. INTRODUCTION 

When a phase difference, A(f, exists between two super¬ 
conductor (S) electrodes separated by a non-superconducting 
material in a Josephson junction, the corresponding charge su¬ 
percurrent is directly controllable via A(f. Motivated by the 
interplay between ferromagnetism and superconductivity, re¬ 
searchers are also interested in the dc Josephson effect in su¬ 
perconducting junctions that contain a central ferromagnet (F) 
region, which in turn can give rise to an additional spin de¬ 
gree of freedom. More specifically, this kind of Josephson 
effect provides a venue for the study of spin currents that can 
be manipulated in cryogenic spintronic systems.Besides 
numerous practical applications^ involving these S FS -based 
Josephson junctions, it is found that novel and interesting phe¬ 
nomena can arise. For example, the realization of a tt state,'^"^ 
where the ground state of the system corresponds to Af = n 
across the junction. Moreover, if the ferromagnet region con¬ 
sists of at least two F layers that each have a uniform magne¬ 
tizations (e.g., a SFFS structure), manipulation of the angle 
between the magnetization vectors can serve to generate long 
range triplet supercurrents^"^"^ in addition to the ordinary sin¬ 
glet ones. Additional control of the magnetic state can also 
occur from the spatially varying spin current within the F lay¬ 
ers of the junction, causing mutual torques to act on their re¬ 
spective magnetic moments. Therefore, 5^5-based junctions 
that contains multiple F layers, present many opportunities for 
controlling the charge and spin currents, and their influence on 
the magnetization in terms of the torque they produce. 

Although interest in the study of F/S multilayer structures 
has recently increased considerably, work on S FS Josephson 
junctions actually started long ago. The Josephson and crit¬ 
ical current oscillations were found to occur as a function of 


the ferromagnet exchange fleld^^’^^ and the thickness of the 
magnet^^. An essential principle behind many of these impor¬ 
tant phenomena is the damped oscillatory nature of the singlet 
Cooper pairs in the ferromagnetic regions, and the associated 
phase shift in the superconducting order parameter. Due to 
the intrinsic exchange fleld in the ferromagnet regions, elec¬ 
trons of Cooper pairs with spin-up (down) decrease (increase) 
their kinetic energy and the Cooper pairs acquire a nonzero 
center-of-mass momentum. It further leads to an oscillatory 
order parameter in the F regions^^. Owing to this oscilla¬ 
tory nature, not only the Josephson critical current but also 
the superconducting critical temperatures oscillate as a func¬ 
tion of exchange fleld and thickness of magnets. Josephson 
junctions of the n type can also be realized by using this prin¬ 
ciple by adjusting either the exchange fleld, the magnet width, 
or both^’^^’^^. The proximity effects between the S and F re¬ 
gions thus give rise to phenomena^^’^^"^^ that subsequently 
play crucial roles in the charge and spin currents that may be 
manipulated in low temperature nanoscale devices, including 
nonvolatile memory elements, where the dissipationless na¬ 
ture of the supercurrent flow offers reduced energy loss and 
Joule heating. 

In equilibrium, singlet Cooper pairs carry no net spin, hence 
any spin current in the system either flows only within the fer- 
romagnets due to their exchange interaction, or it flows by 
means of induced equal-spin triplet correlations, where the 
Cooper pairs have a net spin of m = +1 on the spin quan¬ 
tization axis and they can reside in both S and F regions. 
The generation of long-range triplet proximity effects in su¬ 
perconducting heterostructures with magnetic inhomogene¬ 
ity has been both theoretically predicted and experimentally 
confirmed: by introducing magnetic inhomogeneity, e.g., in¬ 
clusion of an additional magnet with misaligned exchange 
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field, the Hamiltonian no longer commutes with the total 
spin operator and equal-spin triplet correlations can then be 
induced.Due to the imbalance between majority and 
minority spins in a ferromagnet, conventional singlet super¬ 
conducting correlations do not survive for long once inside 
the magnetic region. However, Cooper pairs with electrons 
that carry the same spin are not subject to paramagnetic pair¬ 
breaking and can in principle propagate for large distances in¬ 
side the ferromagnet, limited only by coherence breaking pro¬ 
cesses. Such equal-spin triplet correlations thus play 

an important role in Josephson junctions with inhomogeneous 
ferromagnets.^^’^^ Indeed, it has been reported experimentally 
that with the presence of magnetic inhomogeneity, the Joseph¬ 
son critical current decays much more slowly with increases 
in the F layer thicknesses,^^’^^ as compared to junctions with 
homogeneous magnetization. One of the simplest ways to in¬ 
troduce magnetic inhomogeneity in a Josephson junction, is 
through the insertion of bilayer or trilayer of uniformly mag¬ 
netized ferromagnets. Experimentally, this has the advantages 
of reproducibility and easy manipulation of the relative ex¬ 
change field orientation. These structures also provide direct 
evidence of triplet correlations.Recently, long-range co¬ 
herent transport of triplet pairs was studied in double magnet 
SF 1 F 2 S junctions,^’^^ further demonstrating that it is not al¬ 
ways necessary to have a trilayer^^ ferromagnet structure of 
misaligned ferromagnets to generate equal-spin triplet compo¬ 
nents to the supercurrent. This finding was in consistent with 
the long-range phenomena found in a similar structurewith 
asymmetric widths,and orthogonal exchange fields. The be¬ 
havior of the triplet amplitude is often anticorrelated^^ to that 
of the critical temperature and associated singlet correlations, 
suggesting singlet to triplet conversion. 

The oscillatory and long-range pair correlations also lead 
to new behaviors in the current-phase relations (CPRs) in 
S FS juncitons'^’^’^'^ with a nontrivial magnetic structure. The 
CPR can in general contain not only the first harmonic but 
also higher order harmonics, i.e. I(A(p) ^ I\ sin(A^) -l- 
h sin(2A^). The appearance of additional harmonics in the 
current phase relation has also been discussed in the diffu¬ 
sive and clean regimes for ferromagnetic Josephson junction 
structures.In conventional Josephson junctions with¬ 
out any magnetic interaction, the magnitude of h is much 
smaller than that of Ii . However, in S FS Josephson junc¬ 
tions, near a 0 to tt phase transition, the roles of the first 
harmonic and the second harmonic can be reversed, and the 
CPR can be largely dominated by the second harmonic^^. In 
this regime, both A(f = 0 and A(p = n states can be stable 
or metastable, and they can coexist. The physical origins of 
the second and higher order harmonics are believed to lie in 
the long ranged triplet component of the supercurrent. In this 
regard, within the vicinity of the O-tt transitions, the triplet 
correlations can be tuned accordingly. Since the first har¬ 
monic is suppressed due to the supercurrent fiow reversing 
direction, the higher order harmonics are revealed at the O-tt 
transition point. The influence of interface scattering on the 
higher order harmonics were investigated in the quasiclassi- 
cal clean limit^^ and experimentally detected.^^ The measured 
supercurrent at the O-tt transition point^ was attributed to the 


presence of higher harmonics. Subsequent work with ferro¬ 
magnetic Josephson junctions demonstrated that the higher 
harmonics can naturally arise when varying the location of 
domain walls^^, and also in ballistic double magnetic Joseph¬ 
son junctions, provided that the thicknesses of the magnetic 
layers are unequal^^. Recently, evidence of higher harmonics 
has been observed in Josephson junctions with spin dependent 
tunneling barriers. 

The interaction of the spin current with the magnetization 
in layered ferromagnetic junctions with multiple ferromagnets 
has important consequences for memory technologies. In¬ 
deed, storage of information bits depends on the precise rela¬ 
tive orientation of the magnetizations in two F layers, where 
nonconserved spin currents reflect the mutual torque acting 
on the magnetic moments. The corresponding spin transfer 
torque (STT) can also switch magnetizations when a spin- 
polarized electrical current fiows perpendicular to the layers. 
Spin-transfer torque is known to occur in a very broad vari¬ 
ety of materials, making it an attractive switching mechanism. 
For equilibrium spin currents, governed by spin-polarized An¬ 
dreev bound states'^ ^ , tuning the supercurrent (via A(f ) directly 
infiuences the STT when varying the relative in-plane magne¬ 
tization angle.'^^ The direction of supercurrent flow however is 
not simply related to the direction of the induced torque that 
tends to align the magnetic moments. The triplet correlations 
generated in these types of Josephson junctions (with non- 
collinear relative magnetizations), can also induce spatial vari¬ 
ations in the spin currents responsible for the mutual torques 
acting on the ferromagnets.^^ 

When considering superconducting proximity effects, it is 
important to make sure that the self-consistency condition for 
the pair potential A(x) is fulfilled^ ^ in order to obtain the cor¬ 
rect physical picture. The self-consistency condition is of¬ 
ten neglected in the literature"^^ mainly because it is difficult 
to properly implement in the theoretical studies. As we will 
show in the Sec. II, a source term in the continuity equa¬ 
tion for the charge currents usually arises when a non self- 
consistent superconducting order parameter is used. More im¬ 
portantly, when the self-consistent condition is achieved, the 
free energies of these proximity-coupled systems are properly 
minimized. This concept is crucial especially when study¬ 
ing Josephson junctions, since the superconducting proximity 
effects are the fundamental mechanism behind the nontrivial 
charge and spin currents that fiow within these structures. Fur¬ 
thermore, when the solutions are self-consistent, the charge 
conservation law is satisfied by properly accounting for the 
proximity effects and transport properties for the charge cur- 
rent"^"^. When magnetic inhomogeneity is present, the spin 
density is not conserved and STT can arise and interact with 
the charge dependent quantities. Although there is no con¬ 
tinuity equation for the spin density, since its gradients are 
nonzero, there is still a corresponding, and fundamental, con¬ 
servation law that balances current gradients and the STT^^’"^^. 

Above, we have discussed various physical phenomenon 
associated with the magnetic inhomogeneity such as triplet 
correlations and the generations of STT. Thus understand¬ 
ing the interplay between spin and charge transport with the 
long-range proximity effects is an important topic and consti- 
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FIG. 1. Schematic of the 5'F 1 F 2 F 3 5' Josephson junction considered 
in this paper. A generic configuration is shown, described by the 
in-plane magnetization angles, 0 /, and the out-of-plane angles. Of 
(i = 1,2,3). The ferromagnetic exchange field hf in each region 
is expressed as hi - hi(cos Of, sin Of sin 0 /, sin 6i cos 0 /). 


tutes the main goal of this work. In this work, we therefore 
consider nanoscale S FS Josephson junctions, where the F 
region contains multiple layers: SFFS spin valves consist¬ 
ing of two metallic ferromagnetic layers separated by a non¬ 
magnetic normal metal spacer, and trilayer S FFFS junctions. 
In each scenario, a supercurrent is established via a phase dif¬ 
ference Aif between the S terminals, and different relative 
magnetization orientations are considered. The paper is or¬ 
ganized as follows: We describe our method in Sec II, and 
derive conservation laws for both spin and charge currents. 
Our method is based on solving the microscopic Bogoliubov- 
de Gennes (BdG) equations self-consistently. All important 
physical quantities, e.g., the magnetization, can be extracted 
or constructed from the self-consistent solutions. In Sec III, 
we present a detailed study of the transport properties. We 
conclude with our main findings in Sec IV. 


II. METHODS 

The general method that we use in this paper is that of nu¬ 
merical diagonalization of the self-consistent Bogoliubov-de 
Gennes (BdG) equations. Since many aspects of this method 
have been extensively discussed in previous work^^’^^’^^’"^^ we 
will only include here a brief review of these points, as needed 
to make this paper understandable. We will discuss in more 
detail additional aspects needed for the transport calculations 
described in this work. 

The derivation of the BdG equations for general magneti¬ 
zation configurations begins with the effective BCS Hamilto¬ 
nian, 

= J" d^ry/\r)[F{e-h • cr]i//(r) -f A(r)0^j(r)0^J(r) 

+ A*(r)i/r|(r)i/r|(r)}, ( 1 ) 

where i//(r) = are the usual fermionic operators, 

7-4 = -l/(2m)V^ - Ep + U(r), and cr denote the set of Pauli 


matrices. We describe the magnetism of the F layers by ef¬ 
fective Stoner exchange fields h(r) which in our case have 
components in all (x,y,z) directions (see Fig. 1). The spin in¬ 
dependent scattering potential is denoted by U(r), and A(r) is 
the pair potential. 

To diagonalize the effective Hamiltonian, the field oper¬ 
ators 0^1 and 01 are expanded"^^ by means of a Bogoliubov 
transformation: 

= XI (“«T(^)rn - ’ ( 2 a) 

n 

= X {'^niir'yyn + vl^ir)yl), (2b) 

n 

where Una and Vna are the quasiparticle and quasihole am¬ 
plitudes, which are chosen so that the Hamiltonian is diag¬ 
onalized in terms of the fermionic jn operators. Therefore, 
["H, 7 „] = -enjn and = €nyl. Also, the thermal expec¬ 

tation values involving jn and yl are given by the usual Fermi 
functions The anticommutation relations for 0 and 0 ^ 
yield, 

= (^He - K)\l/^{r) - [h^ - ihy\\l/^{r) + A(r)(/fJ(r), 

(3a) 

[iAx(r),‘K] = (H + - [hx + ihy\\l/^(r) - A(r)^/f|(r). 

(3b) 


It is convenient at this point to simplify to the quasi one¬ 
dimensional geometry (Fig. 1) of our problem. Then, by in¬ 
serting (2) into (3) and using the commutation relations, we 
obtain the general spin-dependent BdG equations for this ge¬ 
ometry. 


' -hx + ihy 

hx - ihy 'Ho + h^ 

0 A* 

A* 0 


0 A '|(M„t(x)^ 

A 0 u„iix) 

-("Ho - h^) -hx - ihy v„|(x) 

-hx + ihy -('Ho + K)) lv„|(x) 


^ lt„l(x) 
■" v„f(x) 


(4) 


where are the quasiparticle energies, x is normal to the lay¬ 
ers, which lie in the y - z plane (see Fig. 1), and A(x) is the 
pair potential, to be found self consistently as explained be¬ 
low. Here the single particle Hamiltonian F(o is written. 


The components of the exchange field h in each of the F lay¬ 
ers take the form: 


hi = h(cos 6i, sin 6i sin 0 /, sin 6i cos 00 , ( 6 ) 

where i denotes one of the magnetic layers. We will assume 
that the magnitude of the exchange field is the same in all 
magnetic layers, and that it vanishes elsewhere. The angles Of 
and 0 / will in general be taken to vary from layer to layer. 
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One obtains in the usual way"^^ the self-consistency condi¬ 
tion for the pair potential, using A(r) = Here 

g is the superconducting coupling constant, which vanishes 
in the non-S layers. In our geometry we find, after using the 
Bogoliubov transformation and making use of the appropri¬ 
ate averages such as (ylyn) = fn, the pair potential can be 
expressed in terms of the quasiparticle amplitudes as an ap¬ 
propriate sum over states: 

A(X) = ImHeJin (7) 

n 

where the prime on the sum indicates that only those states 
that have energies within a “Debye energy”, u;/), are included. 

The problem is then solved iteratively: the potential is ini¬ 
tially taken to be Aq in the first S layer and Aq exp(/A^) in the 
second S layer, where Aq is the initial guess of the magnitude 
for the pair amplitudes. The Hamiltonian is then numerically 
diagonalized and the new pair potential is found via Eq. (7). 
Iteration is continued until convergence. The details of the 
procedure, including the way to ensure, in the Josephson cal¬ 
culations, that the phase difference between the right and left 
ends of the sample remains A^, is explained in Appendix A. 

Equilibrium and transport properties in the structures con¬ 
sidered are strongly influenced by the existence of “odd” 
triplet pairs. The existence of such pairs is allowed by con¬ 
servation laws since, unless all of the F layers have magneti¬ 
zations along the same direction, the total spin of the Cooper 
pairs is not a conserved quantity. Because of the Pauli prin¬ 
ciple, these 5’-wave triplet pairs must have wavefunctions odd 
in frequencyor equivalently^'^’^'^ in time. Within the BdG 
framework, the time formulation is much more convenient. 
Accordingly, we will describe the triplet pair correlations via 
the following amplitude functions, in terms of the field opera¬ 
tors: 

Mr, t) 0)) + 0))], (8a) 

f\{r, t) 0)) “ i'l^iir, 0))]. (8b) 

Taking the quantization axis along the z direction, the triplet 
amplitudes, /o and /i, can be rewritten^^’^'^ in terms of the 
quasiparticle amplitudes: 

/o= (9) 

n 

f, = \i2Y,igV+giyn{t), ( 10 ) 


age of magnetism out of the F layers and into the supercon¬ 
ductors. This can be characterized by the local magnetization 
m(r), 

m(r) = -iiB{r]{r)), (11) 

where r7(r) is the spin density operator, 

r}{r) = j/^^(r)crj/r(r), (12) 

and [iB the Bohr magneton. Eor our quasi-ID geometry, we 
can rewrite the components of m in terms of the quasiparticle 
amplitudes: 

m^ix) = -hb WwntW]/" 

n 

- [v„t(^K|(x) + v„|(x)v;^(x)](l - /„)}. (13) 

my{x) = -ijiB Zl[“ n^{x)ul^{x) - Uni(x)ul^{x)\fn 

n 

+ [v«t(^K|(x) - v„|(x)v*^(x)](l - /„)}. (14) 

nizix) = -llB - \u„i{x)\^\fn 

n 

+ [|v„TWl"-|v„i(x)P](l-/„)). ( 15 ) 

We now turn to the appropriate expressions for the cur¬ 
rents. As stressed in the Sec. I, one needs to carefully estab¬ 
lish proper conservation laws when discussing the transport 
properties of the system^^. We first discuss the charge su¬ 
percurrent. In our geometry the charge current has only one 
component, A which depends on the v coordinate. In the ab¬ 
sence of an external magnetic field, the total charge current, 
Jjc{x) = Jx^{x) is found from the standard quantum 

mechanical expression, Jx = (elm){if/^pxi//). This leads to the 
result: 

JxAx) = ^ 

This expression for the current ensures, together with the self 
consistency condition, that charge conservation is satisfied, 
that is, dJxIdx = 0 in the steady state.It is of course 
convenient numerically to rewrite the expression for the cur¬ 
rent in terms of the calculated quasiparticle amplitudes and 
energies. After inserting the Bogoliubov transformations in 
Eq. (2), we can write the total charge current, as given by 
Eq. (16) summed over spins as: 


where f„(0 = cos(6„0 - /sin(6„0 tanh(6„/2r), and we define 
= UncrV^fj-^- It is sometimes necessary to evaluate the 
triplet amplitudes along a different spin axis. Eor example, 
one may wish to use the direction of the local magnetization 
(defined below) as the axis of quantization. To do so one ro¬ 
tates the quantization axis so that it is aligned with the local 
magnetization direction using the spin rotation matrices dis¬ 
cussed in Appendix B). 

We will consider here spin currents, as well as charge cur¬ 
rents. In our structures spin transport is influenced by the leak- 


Mx) =-Yj (1 - fn)]. (17) 


du 


dvl 


m ^ 

n,(T 


dx 


dx 


One can verify once again the conservation law by taking the 
divergence of the current in Eq. (17) and using the BdG equa¬ 
tions (4), to find:^^’^^ 


dJxjx) 

dx 


2elm I A(x) v„x + v„t ] tanh ( ^ ) 


(18) 
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When the self-consistency condition is satisfied, the right hand 
side vanishes, and charge is properly conserved. If the self- 
consistency condition is not strictly satisfied, the terms on the 
right side act effectively as sources or sinks of current.^ 

We will consider large S contacts with the amplitude and 
phase of the order parameter determined self-consistently ex¬ 
cept near the sample edges (see Appendix A) where sources 
and sinks of charge exist via the implicit external electrodes. 
This gives the necessary charge conservation condition in the 
region of interest. We emphasize here that with self-consistent 
solutions, we are able to correctly determine the effect of 
triplet correlations on both the charge and spin transport. 

The extension of the above considerations to spin transport 
is relatively straightforward.As in the case of the charge 
density, the Heisenberg picture is utilized to determine the 
time evolution of the spin density, r](r,t), 

= (19) 

ot 

where r] is given in Eq. (12). The associated continuity equa¬ 
tion now reads. 


^{r]{r,t))+ — =r+Js, (20) 

Ot ox 

where S is the spin current which in our geometry is a vector 
(in general it is a tensor). The spin-transfer torque, r, is given 
by: 


current in terms of the quasiparticle amplitudes to obtain: 






(27) 


In the case of F layers with uniform magnetization, there is no 
net spin current. The introduction of an inhomogeneous mag¬ 
netization texture however results in a net spin current imbal¬ 
ance that is finite^ ^ even in the absence of a Josephson current. 
This will be discussed in greater detail below. 

To compute the spin transfer torque, it is useful to ex¬ 
press it in terms of the quasiparticle amplitudes. A conve¬ 
nient approach involves directly taking the expectation values 
ofEq. (21): 


T = 2{il/\r)ail/{r)) xh = - m x h, (28) 

Fb 


T = • a, cr]i//ir)} = 2(t//\r)[cr X h]i//(r)}. (21) 

The 2Ts term has components, 

Js. = 21m {A(.Al(r-)«Al(r-) - >lf\{r)>l,\{r))], (22) 

Jsy = 2Re {A(-Al(r)-Al(r) + -A|(r)-A|(r’))} • (23) 

Eor the 5'-wave superconductors considered in this paper, we 
have J's = 0 by virtue of the Pauli principle, since only equal 
time correlations are involved. Thus in the absence spin trans¬ 
fer torque, we have dr]/dt-\-dS/dx = 0. However, in general^ ^ 
spin transfer torque is present and, in the steady state, the 
derivatives of S with respect to v do not vanish. 

The expression for the spin-current, S, is found from taking 
the commutator in Eq. (19) and using Eq. (1): 


where we have used Eq. (11). The magnetization compo¬ 
nents are given in Eqs. (13)-(15). Since the exchange field 
h is prescribed, it is the self consistently calculated magne¬ 
tization that determines the torque acting on the ferromagnet 
layers. Equivalently, one can use the continuity equation in 
the steady state to determine the torque transfer by evaluating 
the derivative of the spin current as a function of position: 


It is however safer to evaluate both sides of Eq. (29) inde¬ 
pendently and use this equation as a consistency check. We 
have performed extensive numerical checks using this proce¬ 
dure. In most of the results presented we have calculated the 
torques using Eq. (28), thus avoiding the numerical derivatives 
that arise when using the right side of Eq. (29). Additional 
physical insight can be gained by integrating Eq. (29) over a 
particular region, e.g., Fi : 


S = - 



dx 


dif/\r) 

dx 


(Ttj/ir)), 


(24) 


SAb) - SAa) = 


I dXTx — 

JFi 


(30) 


where we recall that the vector S represent spin current flow¬ 
ing along the v direction for our quasi-one-dimensional sys¬ 
tems. We can now expand each spin component of the spin 


which means that the change in spin current through Fi (from 
X = aio X = h) equivalent to the net torque acting within 
those boundaries. 




III. RESULTS 
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The results of our systematic investigations are presented 
below in terms of convenient dimensionless quantities. Our 
choices are as follows: all length scales, including the po¬ 
sition X = kfX, and widths Dft = kpcipi (i = 1,2,3) are 
normalized by the Fermi wavevector, kp. For the supercon¬ 
ducting correlation length ^ we choose the value kp^ = 100, 
and the computational region occupied by the S electrodes 
corresponds to a width of 8^ (see Appendix B for numeri¬ 
cal details). All temperatures are measured in units of Tco, 
the transition temperature of bulk S material, and we consider 
the low temperature regime, T/Tco = 0.01. Energy scales 
are normalized by the Fermi energy, sp, including the Stoner 
field interaction h and the energy cutoff, ojp,, used in the self- 
consistency condition, Eq. (7). The latter is set at 0.04: results 
are independent of this cutoff choice. As mentioned above, 
the strength of the magnetic exchange fields, h is taken to be 
the same in both magnets: we set its dimensionless value to a 
representative h = 0.1. We vary the orientation angles of the 
magnetic exchange field in each of the F regions, depending 
on the quantity being studied. The magnetization is normal¬ 
ized by where rie is the electron density, rie = kp/iSn^). 
The normalization tq for the torque follows from the normal¬ 
izations for h and m and Eq. (28). When presenting results 
for the currents, we normalize the charge current densities Jx 
by Jo, where Jo = eUgVp, and vp = kplm is the Eermi ve¬ 
locity. All three components of the spin current S are nor¬ 
malized similarly, by the quantity ^o, where involves the 
normalization of m and a factor of UeVp. The interface scat¬ 
tering U (v) is represented by delta functions of strength H at 
all the interfaces. The corresponding dimensionless param¬ 
eter is Hb = H/vp. The self-consistency of the pair poten¬ 
tial that characterizes an accurate representation of the Cooper 
pair correlations throughout the system, is associated with the 
proximity effects and depends to varying degrees on the pa¬ 
rameters outlined above. In some cases the dependence is 
rather obvious: for instance, large Hb results in weaker prox¬ 
imity effects. In other cases it is more intricate and will be 
analyzed more carefully. 


A. Current-Phase relations 

We begin by showing our results for the self-consistent 
current phase relations in a simple double layer ferromag- 
net Josephson junction. At the interfaces between the F and 
S regions, quasiparticles undergo Andreev and conventional 
refiections.^^"^^ The superposition of these waves in the F re¬ 
gions results in subgap bound states that contribute, together 
with the continuum states, to the total current flow. In Eig. 2, 
we show the supercurrent as a function of phase difference 
(current phase relation CPR) for two ferromagnets of unequal 
width: Dpi = 10 and Dp 2 = 100. This asymmetric choice of 
widths helps ensure^^’^^ that equal-spin triplet correlations are 
generated in the system. The angular parameters in this figure 
are fixed at 6 i = 90°, and 62 = 90°, corresponding to in-plane 
magnetization orientations. The Fi layer has its magnetiza- 



FIG. 2. Normalized (see text) Josephson current versus phase differ¬ 
ence, Aip, for a SF 1 F 2 S structure with Dpi = 10, Dpi - 100, and 
h - 0.1. For panels (a) and (b), the legend in (b) labels the relative 
in-plane magnetization orientations: parallel (P), antiparallel (AP), 
or normal (N). Two interface scattering strengths are considered: (a) 
Hb - 0, and (b) 7/^ = 1. In panel (c), the magnetization orienta¬ 
tions are fixed in the normal configuration (0i = 0, 02 = 90°), and 
the interface scattering is varied as Hp - 0,0.2,0.5,0.6,0.7,0.8,1 (in 
descending order of peaks). 


tion aligned in the z direction ( 0 i =0°). The first two panels, 
(a) and (b), display three different relative in-plane magneti¬ 
zation configurations in the F 2 layers: parallel (P) (02 = 0°), 
antiparallel (AP) (02 = 180°), and normal (N)^^ (02 = 90°). 
Two different strengths of the interface scattering parameter 
are considered. In (a) there is no interface scattering {Hp = 0), 
while in (b), a rather high rate of scattering is present, with 
Hb = 1. The CPR for the collinear configurations (P or AP) 
possesses the conventional In periodicity, and the supercur¬ 
rent flows oppositely for the two alignments. When the rela¬ 
tive magnetizations are orthogonal to each another, the CPR 
becomes n periodic as revealed by the sawtooth-like pattern 
in (a) or the more sinusoidal behavior in (b), both of which 
change sign at A^ 90°. This is a consequence of the emer¬ 
gence of equal-spin triplet correlations^^’^^’^^"^^ that are absent 
when the exchange fields in the ferromagnets point along the 
same direction. When strong interface scattering is present. 
Fig. 2(b) shows that the n periodic CPR (N case) is substan¬ 
tially diminished, relative to the P or AP cases. This is because 
the proximity effect is weakened, with a resulting reduction of 
the associated equal-spin triplet correlations. To further exam¬ 
ine the effects that interface scattering has on this 7r-periodic 
supercurrent, we consider in Fig. 2(c), the same S F 1 F 2 S junc¬ 
tion with varying degrees of scattering strengths Hp and with 
the relative in-plane magnetizations fixed and orthogonal to 
each another (02 = 90°). Increasing Hp clearly leads to a 
crossover in the CPR from a sawtooth to sinusoidal form and 
to a marked reduction of the supercurrent flow. As this occurs, 
the phase difference A(f yielding the critical current density 
also declines. 

We next consider, in Fig. 3, the effect on the CPR when 
changing Dpi and Dp 2 in a SF 1 F 2 S structure. Panels (a)- 
(e) label the various widths Dpi considered. We keep the 
total junction length in which supercurrent fiows constant, 
i.e.. Dp = Dpi Dp 2 , is fixed to a representative value of 
Dp = 390. In these panels, we consider the same three or- 
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FIG. 4. Normalized Josephson current versus A^ for a SF 1 NF 2 S 
structure. The legend labels the N spacer widths, D^. The relative 
in-plane magnetization angle between the two F layers is 90°. For 
panels (a) and (b), the ferromagnetic layers Fi and F 2 have widths 
Dfi = 10, and Dpi - 100. In panel (a) Hp = 0.5, in panel (b) Hp - 
0 .8, and in panel (c), Hp - 0, with ferromagnet widths Dpi - 10, 
and Dp 2 - 380. In panel (c) the additional dotted curve corresponds 
to Dat = 5, illustrating the sensitivity of the current phase relation to 

Dpi. 


FIG. 3. Normalized Josephson current versus phase difference, A^ 
for a SF 1 F 2 S structure with no interface scattering, Hp - 0. In 
panels (a)-(e) various widths of the F regions are considered. Each 
panel is labeled by the width Dpi of the first ferromagnet layer Fi, 
with the constraint that the sum Dpi + Dp 2 - 390. The legend in (b) 
labels the relative magnetization orientation: parallel (P), antiparallel 
(AP), and normal (N). In (f) the magnetization orientation is in the N 
configuration. Dpi and Dp 2 are equal, and their individual widths are 
given in the legend. 

thogonal magnetization orientations as in Fig. 2. One sees 
that, in the N configuration, the structure with the greatest 
geometric asymmetry, as given by the ratio DfilDf 2 , tends 
to have a more pronounced superharmonic CPR relative to 
the P and AP collinear configurations. Upon increasing the 
width Dpi, the jagged sawtooth peaks at Dpi = 10 become 
smoothed and the /r-periodic CPR is reduced substantially. 
Eventually in panel (e), where both F widths are equal, the 
current no longer undergoes a sign change and the additional 
harmonics in the CPR that previously reflected the existence 
of equal-spin triplet correlations are now drastically modifled. 
Remarkably, for this situation, over most of the A(fi range, the 
linear variations of the currents in the orthogonal and antipar¬ 
allel conflgurations tend to overlap. Moreover, it is also evi¬ 
dent that the currents in the N and AP conflgurations that flow 
opposite to the currents in the magnetically uniform P case, 
are substantially greater. The observed reversal of current di¬ 
rection for a given magnetic conflguration when changing the 
F widths is a direct consequence of the damped oscillatory 
behavior of the singlet and triplet correlations, as will be dis¬ 
cussed below (see also Fig. 8). By comparing Fig. 3(a) with 
Fig. 2(a) and noting that the only difference between them is 
the thickness of the second magnet, one flnds that the magni¬ 
tude of the supercurrent for the N case with thicker F 2 does 
not drop signiflcantly, as it does in the P and AP conflgu¬ 
rations. Such behavior is a signature of the long-range na¬ 
ture of the equal-spin triplet correlations. To show further 


the effects of increased width on the supercurrent, we con¬ 
sider in Fig. 3(f) the symmetric geometry conflguration for 
several equal F layer widths Dp. As shown in the legend, a 
broad range of Dp are considered. It is clear from the flgure 
that the equal-spin triplet correlations are strongly impacted, 
and only the 2;r-periodic supercurrent arises. By increasing 
the ferromagnet widths, the rate at which the current changes 
with Acp tends to decline for wider junctions, and the “criti¬ 
cal” phase difference where the current is suddenly reduced 
becomes smaller. In this panel the ferromagnetic region is 
no longer constrained to have the same total width, hence in¬ 
creasing Dp reduces the overall current flow. 

After these examples of two-magnet Josephson junctions, 
we now consider a trilayer junction, where a nonmagnetic 
normal metal “spacer” separates the two ferromagnets. Such 
spacers are often needed experimentally when it is wished^^ 
to rotate the magnetization in one magnet only. To focus on 
the case where the CPRs have important additional harmon¬ 
ics, we keep the in-plane mutual magnetizations orthogonal 
(with 02 = 90°, and 62 = 90°). In Figs. 4(a) and (b), the 
ferromagnet widths are set aX Dpi = 10 and Dp 2 = 100, 
while the width of the normal metal spacer, Dj^, varies from 
zero to 500 (corresponding to f) < < 5) as indicated in 

the legend. The interface scattering parameter has the value 
= 0.5 for the curves shown in (a), whereas in (b) we have 
Hb = 0.8. In either case, the effect of increasing the spacer 
width is to cause an overall reduction in the supercurrent flow. 
For large interface scattering, the CPR becomes less sensi¬ 
tive to Dj^, as seen by comparing (a) and (b). A wider junc¬ 
tion with Hb = 0 is shown in (c), with ferromagnet widths 
Dpi = 10 and Dp 2 = 380. The details of the usual sawtooth 
CPR reveals that in these cases, the supercurrent flow can be 
quite sensitive to the spacer width, as it abruptly reverses di¬ 
rection at considerably different phase differences, for incre¬ 
mental changes in D^. In this panel, results for an additional 
small value of spacer thickness, D^ = 5 are also shown. Since 
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FIG. 5. Normalized Josephson current versus for a SF 1 F 2 F 2 S 
structure. The legend labels the three relative magnetization di¬ 
rections between the F layers. They correspond in order to the 
configurations labeled N P and AP in the previous figures for Fi 
and F 3 , while the intermediate layer F 2 has out-of-plane magneti¬ 
zation. The ferromagnetic layer F 3 has a fixed width correspond¬ 
ing to Df 2 = 380. The adjacent F layers have widths (panel (a)) 
Dfi = Dp 2 - 10 , and (panel (b)) Dpi - Dp 2 - 20 . 


the TT-periodic CPR is closely related to the generation of the 
equal-spin triplet correlation, one can infer from Fig. 4 that 
the introduction of an additional non-magnet metallic layer 
can quantitatively changes the transport properties. 

To examine the effects of increased magnetic inhomogene¬ 
ity, we show in Fig. 5 results for a pentalayer SFiF 2 F^S junc¬ 
tion. This structure is complementary to that studied in Fig. 3, 
the main difference being an additional ferromagnet layer F 2 
with an out-of-plane magnetic exchange field oriented in the 
X direction (corresponding to 62 = 0° in Fig. 1). The relative 
magnetic orientations are labeled in the legend by the direc¬ 
tions of the axes along with the magnetizations are aligned 
in each F layer. For example, z/x/y denotes a sample in 
which Fi and F 3 are normal to each other (the configuration 
labeled N in the previous figures) with an additional out-of 
plane magnetization in F 2 . The width of the F 2 layer is iden¬ 
tical to that of the F\ layer: Dpi = Dp2 = 10 in (a), and 
Dpi = Dp 2 = 20 in (b). The thicker ferromagnet F 3 has 
width Dp 3 = 380. The 7 r-periodic CPR that arises in double 
magnet S F1F2S junctions with orthogonal in-plane magne¬ 
tizations remains relatively unchanged by the addition of the 
additional out-of-plane intermediate ferromagnet. However, it 
was shown in Figs. 3(a)-(e) that the collinear P or AP magnetic 
states with still In periodicities behaved in an approximately 
piecewise linear fashion and that the current maintained its 
direction when varying A(p. Now, on the other hand, the inser¬ 
tion of an additional F layer between the collinear ferromag- 
nets, makes it possible for equal-spin triplet correlations to 
be generated and the supercurrent becomes drastically modi¬ 
fied. For phase differences in the vicinity of 0 or tt, Figs. 5(a) 
and (b) show that the current is approximately linear in the 
phase difference, but there is a broad intermediate range of 
A(f, where the supercurrent flow is relatively uniform. More¬ 
over, for either the P or AP configuration, varying the phase 
can result in the Josephson current switching direction. These 
trends are the same for each of the cases presented in (a) or 
(b). 



FIG. 6 . Normalized Josephson current versus the in-plane relative 
magnetization angle = 0i + 90° (see text), for slSFiF 2 S junction 
(panel (a)), and for rni SFiNF^S junction with Hp = 0.8 (panel (b)). 
The legends show the geometric and material parameters that are 
varied. In (a) and (b) A(p = 100°, Dpi = 10, and Dp2 - 100. In 
(c) an SF 1 F 2 S junction is considered with Hp - 0 and A(p - 60°. 
The outer ferromagnet layer has Dp 2 - 380, while three different 
Dpi values are considered (see legend). 


B. Magnetic orientation and CPR 

Having discussed the current phase relations for a few dif¬ 
ferent ferromagnetic Josephson junction configurations, we 
now will study in more detail the effect of varying magne¬ 
tization orientations on the CPR. We will set the macroscopic 
phase difference to a prescribed value and study the super¬ 
current response within the junctions for a range of relative 
magnetization orientations. 

Beginning with a basic S F1F2S Josephson junction whose 
phase difference is set at A(f = 100°, we consider in Fig. 6 (a), 
the normalized supercurrent density as a function of the in¬ 
plane magnetization angle (p\ {61 = 90°), where we introduce 
0 ' = 01 -f 90°, with 01 being the angle shown in Fig. 1. In 
terms of 0 ^, the P and AP states then correspond to 0 j = 0°, 
and 0 j = 180° respectively (since the magnetization in F 2 is 
fixed along y). Four interface scattering strengths are consid¬ 
ered in panel (a), as indicated in its legend. In all cases, by 
tuning the relative alignment angles, supercurrent switching 
occurs when the mutual magnetization orientations are ap¬ 
proximately orthogonal. As expected, the supercurrent flow 
is greatest for transparent interfaces {Hp = 0 ), and decreases 
with increasing Hp, as the sensitivity to 0 ^ becomes weaker. 
It is notable that the maximum current flow occurs at differ¬ 
ent 0 j values, depending on the interface scattering strength. 
In Fig. 6 (b), an additional N layer, of variable strength as in¬ 
dicated, is inserted between the two ferromagnets and the in¬ 
terface scattering parameter is set io Hp = 0.8. The solid 
= 0 curve corresponds to the = 0.8 curve in panel (a). 
The flgure shows that increasing the N layer thickness tends 
to generally dampen the current through the junction. The 
current flow peaks when the two ferromagnets are aligned in 
the P state and it vanishes altogether near the N conflgura- 
tion, before reversing direction as the relative magnetizations 
approach the AP state. In Fig. 6 (c) we investigate the super¬ 
current flow in a wide S F1F2S junction, with Dp2 = 500, and 
three different Fi widths (see legend). The macroscopic phase 
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FIG. 7. Normalized Josephson current for a SF 1 F 2 F 2 S junction 
versus (panel (a)) the out-of-plane angle 62 (see Fig. 1), and (panel 
(b)) the in-plane angle A set phase difference Acp = 100° is as¬ 
sumed for all cases. In (a) the ferromagnets, Fi and F 2 have widths 
Dfi = Df 2 = 10 , and the legend lists the three F 3 widths considered. 
The exchange field points along the j direction in F 3 and the z direc¬ 
tion in Fi. The right vertical axis is for the Z )^3 = 100 case only. For 
the Df 2 = 380 case, there is no interface scattering (Hb = 0), while 
the remaining Dpi, cases have Hb = 0.8. In (b) 61 = 90°, and h is 
directed along y in F 2 , and along z in F 3 . The legend indicates the 
F layer widths of the two trilayer ferromagnet structures considered: 
they have the same total width. Interface scattering in both cases is 
Hb = 0 . 8 . 


difference is set at A(p = 60°. It is seen that for this geome¬ 
try, the current has some contrasting features compared to the 
previous cases involving thinner F layers. For instance, when 
tuning 0 ', now the weakest current flow occurs when in the 
P state, and the maximum occurs in the orthogonal configura¬ 
tion. For all F width considered, the current undergoes rapid 
changes in the vicinity of the middle value between N and P 
or N and AP configurations, with iht Dpi = 10 case changing 
the most, and then remains nearly constant at angles near the 
P and AP configurations. Since the prescribed phase differ¬ 
ence A 0 = 60° is near where the maximum currents occur in 
TT-periodic Josephson junctions, it is evident that the results in 
Fig. 6 (c) are a direct consequence of the induced equal-spin 
triplet correlations. We will see below in Sec. IIIC that this 
correlates with the triplet generation behavior. 

We next move to the study of the magnetization orienta¬ 
tion role in the supercurrent for more complicated S F 1 F 2 F 3 S 
junctions. As an example, we consider a scenario where the 
Fi and F 3 layers have magnetizations pinned along the z and y 
directions respectively, that is, a relative N configuration, but 
the magnetization in the central F layer rotates on the xz plane 
(02 = 0) from along the x-axis to the z-axis. In Fig. 7(a), we 
show the normalized current density as a function of 62 for few 
different F^ layer widths, Dp^. We set Dpi = Dp 2 = 10, and 
consider values of Dp^ that lead to both symmetric and asym¬ 
metric structures. In each case, A^ = 100°, and interface scat¬ 
tering is present with Hp = 0 . 8 , except for the widest junction 
with Dp 2 > = 380 where the interfaces are transparent {Hp = 0). 
In the case Dpi = Dp 2 = Dpi, = 10, where all F thicknesses 
are identical, the current flow is approximately antisymmet¬ 
ric around 62 = 180°. The current flow is suppressed for the 
asymmetric Dp^ = 100 situation by decoherence effects aris¬ 
ing from the larger width. While Jx for the simpler S F 1 F 2 S 



FIG. 8 . Normalized Josephson current for a SF 1 F 2 S structure as 
a function of F layer width Dpi. The other ferromagnet’s width is 
fixed Bi Dp2 = 100. The exchange field is along z in Fi, and along 
y in F 2 . The legend indicates two normalized h values: The dashed 
curve corresponds to h = 0.04, while the solid curve is for h = 0. 1 
The interface scattering parameter is set to Hp = 0.8, and a phase 
difference of Aip = 100° is maintained across the S electrodes. 


junctions in Fig. 6 was tt periodic in 0 p variations in O 2 for 
trilayer ferromagnetic junctions are in general In periodic, as 
seen in the figure. The asymmetric case also demonstrates a 
more intricate structure as the magnetization angle is swept. 
For the narrower junction with Dp^ = 10, the orientations 
in which the current switches direction are near 62 = 0 ° and 
O 2 = 180°. When Dp^ = 100, the current never changes its 
direction. Although the rotating out-of-plane exchange field 
of F 2 does affect the strength of current with 2 n periodicity. 
For the highly asymmetric Dpi, = 380 case, the current again 
maintains its direction of flow over the full angular range of 62 . 
Also, the current is approximately constant except for orien¬ 
tations when the exchange field in F 2 points along z: O 2 = 90° 
or 62 = 270°. For orientations near 62 = 270°, we find that 
the current is strongly suppressed. In Fig. 7(b), the in-plane 
magnetization in the first ferromagnetic layer is now varied, 
while the other two are kept fixed: for F 2 , it is along y, and for 
F 3 , it is along z, an N configuration. Two types of structures 
are considered, with the total width of the three ferromagnetic 
regions being constant. In the first case, Dpi = Dp 2 = 10, and 
Dpi, = 100, while in the second one all three F layers have 
width, Dpi = Dp 2 = Dpi, = 40. Both cases exhibit similar 
behavior as a function of 0 ^. The current vanishes when Fi is 
antiparallel to F 2 and is highest when the magnetization lies 
nearly in between the those of F 2 and F 3 . Thus, the charge 
supercurrent which flows oppositely in the two structures, can 
be effectively switched on or off by manipulating the in-plane 
magnetization angle of the first ferromagnet. 

It is known that in FjS heterostructures, including 
bilayers^^ and Josephson junctions, variations in the magnetic 
exchange field and ferromagnet thickness induce damped os¬ 
cillations in the spatial behavior of the Cooper pair ampli¬ 
tudes, resulting in modulation of physical quantities as a 
function of either h or Dp. The damped oscillations in the 
clean limit have a wavelength that goes as the inverse of 
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FIG. 9. Normalized triplet correlations in an SFiNF^S Josephson 
junction as a function of position. The F layers have widths Dpi = 10 
and Dp 2 = 380, and Hp - 0, corresponding to the parameters used 
in Fig. 4(c). The top set of panels relate to structures with a normal 
metal spacer, N, of width = 10, while the bottom set represent a 
larger Allayer with T)AT = 100. The dashed vertical lines represent the 
interface between the N and F 2 regions. Various phase differences 
A(p are considered (see legend). The magnetization in Fi is along z, 
while it is along y in F 2 . 


the exchange field. To investigate this phenomenon in an 
SF 1 F 2 S junction, we show in Fig. 8, the supercurrent that 
fiows through the junction as a function of the Fi width, Dpi. 
For the adjacent F 2 layer, the width is 0^2 = 100, and with 
interface scattering strength Hp = 0.8. Two exchange fields 
are considered: h = 0.04 and h = 0.1. The period of os¬ 
cillations in each case is seen to be approximately Inkp^r, 
where kp^p = sp/h = 25, and 10, respectively. The current 
for both cases is maximal when the two F regions are equal, 
and then slowly dampens out with increasing Dpi. This decay 
length is inversely proportional to the magnitude of the ex¬ 
change field. Moreover, the charge current in each case peri¬ 
odically changes sign. Since the oscillations in the current as a 
function of Dpi increase with larger exchange fields, the cur¬ 
rent direction can be very sensitive to fabrication tolerances 
for strong magnets. 


C. Induced triplet pairing 

We now discuss the induced triplet pairing correlations in 
ferromagnetic Josephson junctions. The presence of multiple 
misaligned ferromagnets yields both the m = 0 (Eq. (9)) and 
the m = +1 (Eq. (10)) triplet pair amplitudes as permitted by 
conservation laws and the Pauli principle. To gain an over¬ 
all view of the opposite-spin triplet amplitudes, /o, and the 
equal-spin amplitudes, /i, we illustrate in Eig. 9 the spatial 
behavior of these correlations in the N and F regions of an 
SF 1 NF 2 S junction. We focus on the real parts of these com¬ 
plex quantities, keeping in mind that the imaginary compo¬ 


nents obey similar trends. The geometrical parameters in this 
figure are Dpi = 10, Dp 2 = 380, and D^ = 10 (top panels) 
or Dn = 100 (bottom panels). Thus, in all panels the region 
800 < X < 810 is occupied by Fi while F 2 occupies the re¬ 
gion 820 <X < 1200 in the top panels and 910 < X < 1290 
in the bottom panels (c) and (d), where the vertical dotted line 
denotes the N spacer boundary. Hence, different horizontal 
scales are used in each case. The scattering parameter is set to 
Hb = 0 , and each curve corresponds to a different phase dif¬ 
ference, as indicated by the legend. The exchange fields 
are in-plane and normal to each other. Within the Fi region, 
panels (a) and (c) reveal that the magnitude of the /o pair cor¬ 
relations are approximately of the same magnitude, decreas¬ 
ing in the vicinity of the SIF interface located at X = 800. 
The system with the wider normal metal layer, D^ = 100 is 
slightly more sensitive to phase variations. Within the ferro- 
magnet F 2 , the same panels (a) (c) show the oscillatory nature 
of /o, which behaves similarly to the singlet pair amplitude, 
the periodicity arising from the difference in spin-up and spin- 
down wavevectors. Eor the chosen exchange field, the oscil¬ 
lations are limited in Fi due to the confined width. Turning 
now to the equal-spin triplet correlations /i, panels (b) and 
(d) display behavior which contrasts with the /o results. In 
particular, within the narrow Fi region, the fi triplets are neg¬ 
ligibly small, and the /o correlations clearly dominate. In the 
N region, the /o correlations nearly vanish, while the equal- 
spin triplets peak near the Fi/N interface (at X = 810), be¬ 
fore dropping within the normal metal (see Eig. 9(d)). Einally, 
within F 2 , the triplets fi assume a slow, long-range variation 
compared to the damped oscillatory behavior of the /o curves. 

In addition to investigating the spatial behavior of the triplet 
amplitudes, it is instructive to also examine the spatially av¬ 
eraged triplet and singlet correlations as functions of the rel¬ 
evant system parameters. Eor example, by tuning the relative 
magnetization angle, (varying 0 ' at fixed 02 = 90°), impor¬ 
tant overall features can be revealed. The top panels (a) (b) 
(c) of Eig. 10 show the 0^-dependence of the magnitudes of 
the triplet and singlet amplitudes averaged over the F 2 region 
for the S F 1 NF 2 S structure studied initially in Eig. 6 (b). Eour 
representative N layer widths are considered as indicated in 
the legend. The proximity effects and hence coupling of the 
two ferromagnets diminish with increasing D^, and therefore 
the pair correlations become less sensitive to variations in 0 ', 
as observed for the largest D^ = 500 case. Other than the 
diminished magnitudes, the overall trends and behavior how¬ 
ever do not depend strongly on the presence of the normal 
metal spacer. This may be important in experiment design, 
where spacers are often needed. The opposite-spin triplet cor¬ 
relations, /o and the singlet pair amplitude (/s = A/g) in (a) 
and (c) behave in rather similar ways, but /s is more symmet¬ 
ric about the orthogonal direction 0 ' = 90°. When the rela¬ 
tive magnetization orientation varies in inhomogeneous S/F 
systems, the process of singlet-triplet conversion plays a role 
in the transport and thermodynamic properties of such sys¬ 
tems. It is apparent from panels (a) and (b) that the orienta¬ 
tion 0 ' that leads to a minimum in /o and /s, corresponds to 
that where fi is largest. These occurrences arise when the ex¬ 
change fields in Fi and F 2 are nearly orthogonal. This angle 
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FIG. 10. Top row, panels (a)-(c): Normalized singlet and triplet cor¬ 
relations versus in-plane relative magnetization angle (p[ = 0i + 90° 
for SLSF 1 NF 2 S structure. The magnitudes of these pair correlations 
are averaged over the F 2 region. The geometrical and material pa¬ 
rameters are the same as in Fig. 6(b). The bottom panels (d) and (e) 
correspond to the local spatial behavior of the triplet correlations for 
the Dn = 100 case studied in the top panels. The two F regions 
are plotted in separate frames to discern the triplet correlations in the 
narrow Fi region. The dashed vertical lines represent the interface 
between the N and F 2 regions. Several values of (p[ are considered, 
as shown in the legend. 


for /o slightly shifts with variations in D^. 

In the bottom panels (d) and (e) of Fig. 10 we display the 
spatial dependence of the real components of the triplet cor¬ 
relations throughout each of the three junction regions dis¬ 
cussed in the top panels. Results for five different relative 
magnetic orientations are presented, (see legend) including 
the P ( 0 ; = 0°), AP ( 0 ; = 180°), and N (0; = 90°) con¬ 
figurations. Considering first Fig. 10(d) in the Fi region, we 
see behavior similar to that found for the wider junction case 
in Fig. 9, including a relatively weak dependence on the orien¬ 
tation angle 0j (as opposed to the phase difference A(f). The 
central N region is most affected by variations in 0j: the real 
part of /o changes sign when 0j is swept from the P to AP 
state, and nearly vanishes altogether at (p[ ^ 90°. In F 2 , 3. 
series of oscillations emerge with a periodicity similar to that 
found in Fig. 9, since an exchange field of the same strength 
was used. In Fig. 10(e), the equal-spin /i amplitudes exhibit 
considerably different behavior. First, only three of the con¬ 
sidered 0j yield nonzero results, since the P and AP configu¬ 
rations cannot generate equal-spin triplet correlations. Inter¬ 
estingly, within the Fi region fi does not exhibit a slow decay, 
but rather oscillates with a period that is much shorter than the 
opposite spin pairs governed by the difference in spin-up and 
spin-down wave vectors. We find that within the normal metal 


FIG. 11. Normalized triplet |/o,avgL l/i,avgL and singlet |/3,avgl am¬ 
plitudes, averaged over the F regions indicated on the overhead leg¬ 
ends, plotted as functions of 62 and (p[ The geometrical and material 
parameters corresponds to the Dp 2 - 100 cases presented in Fig. 7. 


layer, there is nearly a complete absence of equal-spin correla¬ 
tions, this is accompanied by the appearance of opposite-spin 
correlations /o (see panel (d)). The f\ amplitudes are largest 
in the F2 layer, for the relative orientation of (p[ = 90°, in 
agreement with the averaged results in Fig. 10(b). For the 
relative magnetization angles of 0j = 30°, and 0j = 150°, 
the f\ amplitudes are identical due to the symmetry about 
0 ' = 90°. One can now correlate the features Fig. 10 and 11 
to Fig. 4. From Fig. 4, we learn that the N magnetic configu¬ 
ration often leads to the appearance of 7r-Josephson junctions. 
Here we are able to give concrete proof that in the N cases 
the equal-spin triplet correlations are maximized and are in¬ 
sensitive to A^. Therefore, the CPRs for different magnetic 
configurations are essentially characterized by their detailed 
singlet/triplet nature. 

We next examine, in Fig. 11, the behavior of the aver¬ 
aged triplet and singlet amplitudes, as the magnetic orien¬ 
tation angles, 62 (top panels), and 0j (bottom panels) are 
changed, in more complicated SF 1 F 2 F 3 S Josephson junc¬ 
tions. This study is therefore complementary to the results 
shown in Fig. 7 involving the charge supercurrent. The ge¬ 
ometric parameters are Dpi = Dp 2 = 10, and = 100. 
The region in which the pair correlations are averaged over 
is specified in the top legends. In the top row of Fig. 11 we 
present results for magnetization orientations 62 sweeping the 
entire angular range from 0° to 360°, while the magnetizations 
are aligned along z in F\ and along y in F 3 . Therefore when 
62 = 0° or 62 = 180°, all three ferromagnets have mutually 
orthogonal magnetizations, corresponding to a high degree 
of magnetic inhomogeneity. Under these circumstances, one 
can expect that the equal-spin triplet correlations fi should 
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be, on the average, at their highest values, while the opposite 
spin correlations should be weakest. Indeed, in the regions 
7^3 and F 2 , the opposite-spin singlet /s, and triplet /o corre¬ 
lations possess minima near these angles. In contrast to the 
spatially averaged /i amplitudes peak, in F 3 and F 2 , at those 
orientations. Although the general trends are usually the same 
for all, F layers, proximity geometrical effects can result in 
self-consistent triplet correlations with more intricate nontriv¬ 
ial structure, and this is the case with the averages over the 
Fi region. In the bottom set of panels of Fig. 11, we consider 
in-plane magnetization rotations of the Fi layer. The other 
ferromagnets F 2 and F 3 have their magnetizations fixed in the 
y and z directions respectively. For this situation, the opposite- 
spin /o amplitudes in Fi are seen to be 2 n periodic, peaking at 
0 ' 225°. In F 2 these triplet amplitudes are seen to be largest 

when the relative orientations between F\ and F 2 are either P 
( 0 ' = 0°,360°) or AP ( 0 ' = 180°). This is consistent with 
the behavior of the triplet amplitudes found in double magnet 
spin valve systems^ ^ 


D. Spin Transport 

Having established the salient features of supercurrent 
charge transport and pair correlations in a variety of ferromag¬ 
netic Josephson junctions, we now explore the spin degree of 
freedom and determine the crucial spin currents and the asso¬ 
ciated spin transfer torques. The current that is generated from 
the macroscopic phase differences between the S electrodes 
can become spin-polarized^when entering one of the fer¬ 
romagnetic regions. This spin current can then interact with 
the other ferromagnets and be modified by the local magneti¬ 
zations due to the spin-exchange interaction, via the existence 
of spin transfer torques. The conservation law associated with 
process is described by Eqs. (28) and (29). It is important 
not only to understand the behavior of the spin-polarized cur¬ 
rents in ferromagnetic Josephson junctions, but also the vari¬ 
ous ways in which to manipulate them for practical spintronic 
applications. We therefore investigate from a microscopic and 
self-consistent perspective, the equilibrium spin currents and 
associated torques throughout the entire junction regions as 
functions of position, phase difference, and magnetization ori¬ 
entation angles. 

In Fig. 12(a), we consider the spatial dependence of the spin 
current in a SF 1 NF 2 S junction. The geometrical parameters 
used in this plot are the same as in Fig. 4(a), with = 100. 
Our geometry ensures that the generally tensorial spin current 
is reduced to a vector in spin space, representing a spin vector 
current flowing in the spatial v direction and having in general 
three components in spin space. We display the spatial depen¬ 
dence of the V spin component, Sx, (normalized as previously 
discussed). Because the exchange interaction (and hence the 
torque) vanishes in the N and S regions, only the F regions of 
the junction can have a spatially varying spin current: in the N 
and S regions the spin current must be spatially invariant. Un¬ 
der our constant-phase and zero voltage boundary conditions 
the outer 5 ’-wave superconducting regions do not^^ support a 
spin current, and hence S vanishes there. The central nonmag- 


xlO^ 




FIG. 12. Panel (a): Normalized x component of the spin current 
in an Josephson junction vs. position. The geometric 

widths are Dpi = 10, Df 2 = 100, and = 100. There is moderate 
interface scattering with Hb = 0.5 (see e.g.. Fig. 4(a)). The dashed 
vertical lines mark the interfaces. Panel (b) The total torque r^c^tot 
acting within the Fi region as a function of Aip. Several normal metal 
widths are considered, see legend. The = 0 case has been shifted 
downwards by 4 x 10"^ for comparison purposes. In panel (c) the 
^-component of the local torque is shown vs. position with the same 
phase differences used in (a). The inset is a magnification of the 
torque within the narrow Fi region (800 < X < 810). 


netic N layer, however, couples the two ferromagnets via a 
constant spin current, which is related to the net torque acting 
within the F regions (see Eq. (30)). The spin current oscillates 
in the F 2 region. The amplitude of these oscillations depends 
on A(f, while the period does not: the points in F 2 where Sx 
vanish are independent of Aip. If a ferromagnet is very thin, 
as occurs for F\, the spin currents vary nearly linearly with 
X, which can be viewed as a small segment of a sinusoidal 
function. To present an overall view of how the change in 
spin current and its associated torque vary as the phase varies. 
We do this in terms of the total torque, defined as the integral 
of the local torque, normalized as discussed above, over di¬ 
mensionless distance. We plot in Fig. 13(b) the total torque, 
U,tot 5 within Fi. for a few values of the interlayer N spacer 
thickness. In all cases, Tx,tot is relatively uniform until a sharp 
crossover near A(f = 90°, where the net torque changes sign, 
coinciding with the point of supercurrent reversal (Fig. 4(a)). 
Interestingly, only when the normal metal insert is present, 
does Tx,tot reverse direction. Panel (c) illustrates the local v- 
component of the torque, Tj^, as a function of position through¬ 
out the entire junction region. The inset is the same quantity. 
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FIG. 13. Panel (a): Normalized spin current in m SF 1 NF 2 S Joseph- 
son junction as a function of position. The geometric widths are 
Dfi = 10, Df 2 = 100, and = 100. The interface scattering 
strength is set to = 0.8. These values are the same as in Fig. 6 (b). 
The dashed vertical lines mark the interfaces. Panel (b) illustrates 
the total torque Tx,tot within the Fi region as a function of the rela¬ 
tive in-plane magnetization angle, 0 j. Several normal metal widths 
are considered as depicted in the legend. The = 0 results have 
been multiplied by 1/5 for comparison purposes. In panel (c) the 
v-component of the torque is shown as a function of position for the 
same case considered in (a). 


but plotted only over the narrow Fi region. To correlate with 
(a), the = 100 case is considered here. Each curve repre¬ 
sents a different phase difference as shown in the legend for 
(a). For in-plane exchange field interactions, no other compo¬ 
nent of the torque can exist in equilibrium when spin currents 
do not enter or leave the superconducting electrodes.Thus, 
the net torque for the entire system must vanish, requiring Tx^ot 
for each of the two F regions to be opposite in sign, despite 
the dissimilar spatial behavior as exhibited in (c). Compar¬ 
ing (a) and (c), it is also evident that within the oscillatory F 2 
region, Tx and Sx behave similarly, but are out of phase by ap¬ 
proximately 90°, in agreement with Eq. (29). Since the spin 
current was shown to be uniform in the normal metal region, 
the torque is seen to vanish there, as it should be in regions 
where the magnetic exchange interaction is absent. 

We now turn to studying how the spin currents and associ¬ 
ated torques change when varying the relative exchange field 
directions between Fi and F2 in an S F1NF2S junction. A su¬ 
percurrent is generated in the structure by maintaining a phase 
difference which we take to be A(p = 100° in Fig. 13. We ro¬ 
tate in-plane magnetization Fi, while keep that in F 2 fixed 


- F _ F __ F 




A 


FIG. 14. Top panels: total torque within each of the ferromagnet 
regions (see overhead legend) in an junction as the angle 

62 varies. The system parameters are those used in Fig. 7(a). The 
sum T/,tot ii - -T, y, x) over all three ferromagnetic regions vanishes 
for each component. In the bottom panel, the spatial behavior of the 
normalized v-component of the spin current is shown throughout the 
system for a few select magnetization orientations 62 (see legend). 
The inset is a magnification of the region centered around the E 2 /F 3 
interface located at X = 820. Vertical dashed lines in the main plot 
mark interface locations. 


along the y (as in Fig. 6(b)). Control of the free-layer mag¬ 
netization can be achieved experimentally via external mag¬ 
netic fields^^, or spin-torque switching^^’^^. In Fig. 13(a), the 
X component of the local spin current, S xA^ shown through¬ 
out the junction as a function of position X, for four values 
of the (/)[ angle (we have 61 = O 2 = 90°). The spin current 
is again seen to be a nonconserved quantity within the ferro- 
magnets, reflecting the existence of a spin-transfer torque. In 
the nonmagnetic normal metal connecting the two F regions, 
the current is constant, and its value as (f)[ is varied Sx cycles 
from positive to negative. To explore this further, we exam¬ 
ine the total change in spin current across Fi, as 0' sweeps 
from the P to AP state. This change is related via Eq. (30) 
to the integrated torque in this region. Hence, in Figure 13(b) 
we plot Tx,tot vs (p[ for a wide range of D^. When the normal 
metal is absent (D^ = 0), the magnitude of the total torque 
reaches its peak around 0^ = 90°, indicating that this compo¬ 
nent of the torque, which tends to align the magnetic moments 
of the two F layers is largest when they are mutually orthog¬ 
onal. This makes sense physically. The presence of even a 
thin normal metal spacer causes Tx,tot to become much smaller 
(the Dn = 0 results are plotted after dividing them by five) 
and nearly tt- symmetric, so that now the orthogonal magnetic 
configuration produces negligible net torque within the F lay¬ 
ers. Increasing reduces the ferromagnetic coupling and 
hence reduces the magnitude of the mutual torques, although 
the TT-periodicity is retained. Finally in Figure 13(c) we plot 
the local value of the torque, and find its spatial behavior to be 
consistent with that of 5 as given in Eq. (29). 
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Finally, we consider the SFiF 2 Fi,S system, studied previ¬ 
ously in Fig. 7(a), with Dpi = Dpi = 10, and Dp^ = 100. A 
phase difference of A(f = 100° maintains a constant current 
throughout the junction, and there is moderate interface scat¬ 
tering, with Hb = 0.8. The magnetization in Fi is along z, 
and in F 3 , it is along y. The central ferromagnet, F 2 , has a 
magnetization vector that is rotated the vz plane (see Fig. 1), 
so that for ^ = 0°, it is oriented along x, and for 6 = 90°, it is 
aligned along z. For these more complex magnetic configura¬ 
tions, where one of the F layers possesses an out-of-plane ex¬ 
change field, all three spin components of the current S must 
be considered. The top panels in Fig. 14 depict the compo¬ 
nents of the total torque, T;;tot (i = z) for each ferromagnet 
region, identified in the legend above these panels. Since the 
total torque in a given direction equals (see Eq. (30)) the over¬ 
all change in spin current, and, as explained above there is 
no spin current in the S regions at fixed phase, the sum of 
each component T/,tot over all F regions must be zero. This is 
seen in these three panels, where the oscillatory curves exactly 
cancel one another. For each of the three components, we also 
observe that the total torque in either Fi or F 3 nearly vanishes 
over the whole angular range of 62 . This follows from the ex¬ 
pression for the local torque, Eq. (21), which implies that r is 
orthogonal to the exchange field vector h, and the magnetiza¬ 
tion m. For example, considering the leftmost panel the only 
component of the exchange field in Fi is along z, and since 
Tx ~ myhz, a y component of the magnetization is needed in 
Fi to generate a torque along x. However h in the adjacent 
F 2 rotates solely in the vz plane, and thus nty vanishes in Fi 
except in a narrow region near the interface (see Ref. 31) re¬ 
sulting in a very small value for the averaged Tx,tot- 

As mentioned previously the spin current S' is a local quan¬ 
tity, and when it is spatially nonuniform the resulting torque 
infiuences the magnetization configurations. It is therefore in¬ 
sightful to examine also in this case, as we did in Figs. 12 and 
13, the spatial behavior of the spin current. The results are dis¬ 
played in the bottom panel of Fig. 14. For clarity, we present 
only the v-component, Sx, since the other spin components 
behave similarly. A few representative angle orientations, O 2 , 
are considered (see legend). As expected, we see that Sx van¬ 
ishes in the outer S electrodes. Within the larger F 3 region, 
the spin current undergoes regular oscillations which are much 
harder to distinguish in the narrow F 2 and Fi regions. To com¬ 
pare with the previous results, we see from this panel that the 
change in the v-component of the spin current, AS x across the 
F\ boundaries is negligible, in agreement with the results in 
the left top row panel. In F 2 , this component of the current 
is very small near the left interface, but it increases near the 
right edge, so that AS x agrees well with the Tx,tot variations 
in F 2 observed in the left top panel. Also in agreement is the 
enhancement of the spin current S x that is narrowly peaked 
at the F 2 /F 3 interface for O 2 that are near normal to the plane 
{O 2 ~ 90°). Finally, since the interface adjoining F 3 and the 
right S electrode has S x = ^, the change A^j^ in the F 3 region 
is due entirely from the value of the spin current at the F 2 /F 2 
interface, thus resulting in A^j^ that is exactly opposite to that 
in F 2 . 


IV. CONCLUSIONS 


We have presented here an extensive study of the Joseph- 
son currents flowing in generic ballistic structures of the S FS 
type where the F regions contain two or three ferromagnetic 
layers, possibly separated by normal spacers. For the S FFS - 
type spin valves, we study their transport properties by consid¬ 
ering different in-plane relative magnetization angles. When 
the third F layer is present (S FFFS ), we allow the central F 
region to have out-of-plane magnetic orientation while those 
for the two outer F layers are still in-plane. Our self consistent 
formalism ensures^ ^ that the charge conservation law is satis¬ 
fied and that the proper relations that balance the spin transfer 
torques and the gradients of the spin current components hold. 
Results are given for a wide range of values of the geometrical 
and orientation parameters, as well as interfacial scattering. 

We have organized our results in several subsections. We 
first have considered (Sec. Ill A) the current-phase relations 
(CPRs) as a function of geometrical parameters (layer thick¬ 
nesses) at fixed relative angles between the in-plane magne¬ 
tizations: parallel (P), antiparallel (AP), and normal (N). We 
find that in general, larger geometric asymmetry (the aspect 
ratio for the thicknesses of the two outer F layers) leads to 
larger superharmonic (tt periodic) behavior. This is particu¬ 
larly pronounced in the two-magnet case. It is found that the 
strength of interfacial scattering can affect the magnitude of 
the critical current. Next, in Sec. IIIB, we consider the ef¬ 
fect of magnetization misorientation on the CPR. We find that 
these effects are profound. In particular, by sweeping the rel¬ 
ative in-plane angle from P to AP at fixed phase differences 
A(f between two S electrodes, the supercurrent flow first van¬ 
ishes at N configurations, followed by reversal of its direction. 
These can be understood to a very large extent by noting that 
the generation of induced spin triplets (studied in Sec. Ill C) 
is correlated with magnetic inhomogeneity and, via this phe¬ 
nomenon, to the CPR relationships. In Sec. Ill C, we also pre¬ 
sented the local spatial behavior of both the m = 0 and m = ±l 
triplet correlations and carefully quantify each component as 
functions of magnetic orientations and A(p. The results clearly 
demonstrate the existence of the singlet-to-triplet conversion 
in the Josephson junctions. Finally, in Sec. HID we have dis¬ 
cussed spin transport and shown that, due to the interaction 
between the charge current and the magnetizations, both the 
spin current and the spin transfer torque (STT) oscillate in the 
F regions. By varing Acp or the magnetic misalignment angles, 
the phase of their oscillations can change accordingly. In ad¬ 
dition, we have shown how the spin current gradient equals 
the spin transfer torque. 

We hope that our rather comprehensive study of transport 
in these multilayer structures will guide the experimentalist in 
choosing optimal configurations for building devices such as 
low dissipation memory storage units, which are expected to 
rely on the behavior of the Josephson junctions studied here, 
in particular on the behavior of the CPR as orientation angles 
are changed. 
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Appendix A: Numerical Procedure 


Here we discuss some technical aspects of the numerical 
procedure used in calculating the spin and charge currents 
governed by the Andreev bound states. We first expand^^ the 
quasiparticle amplitudes in terms of a complete set: 





^m{kqX)ll/q{kq), 


(Al) 


where we use the shorthand notation = 

{Un^{x),Uni{x),Vn^{x),Vni{x)), and if/q = {Uq^,Uqi,Vq^,Vqi). 

We write the wavevector index as, = qnld, so that 
lS.kq = kq+i - kq = 71 jd. Thus N grid points subdivide 
the system of width d. We take d to be large enough so 
that the results become independent of d. The next step 
involves Fourier transforming the real-space BdG equations 
(Eq. 4), resulting in the following set of coupled equations in 
momentum space: 


Hamiltonian are expanded using Eq. (A5) and Eq. (A3) re¬ 
spectively. We then find the quasiparticle energies and ampli¬ 
tudes by diagonalizing the resultant momentum-space matrix 
(Eq. (A2)). Once the momentum-space wavefunctions and en¬ 
ergies are found, they are transformed back into real-space 
via Eq. (Al). Erom them, a new pair potential A(x) is self- 
consistently determined via Eq. (7) through the entire region 
except for a small region (three coherence lengths from the 
sample edges) where the pair potential is fixed to its bulk ab¬ 
solute value, with phases 0, A^. The newly calculated A(x) 
is then used in the BdG equations and the above process is 
repeated iteratively until convergence is achieved. When de¬ 
termining the current phase relations, A^ is defined as the dif¬ 
ference in phases between the superconductors in the outer¬ 
most self-consistent regions. As self-consistency evolves with 
each iteration, the final A^ often differs slightly from the fixed 
difference A^ that is set in the non self-consistent regions. 
Thus, to have A^ fixed to a prescribed value while varying 
other parameters, e.g., 62 , additional calculations are needed 
with slightly different initial choices for the phase A^. Eol- 
lowing the discussion in the main text, when current is fiow- 
ing through the junction, the self-consistently calculated re¬ 
gions are always been found to possess the necessary spa¬ 
tially constant current. The fixed-phase non-self consistent 
edge regions, provide the physically necessary source or sink 
of current, via the applied electrodes, thus acting as an effec¬ 
tive boundary condition. 
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(A2) 


Here we have defined = (fii^-, U 2 o-, • • •, unct), Vo- = 

(vicr, V 2 cr,..., va^o-), and the matrix elements. 


2 

Ho(q,q') = - 
d 

2 

Kq.q') = 

d 

2 

hi{q,q') = - 

d 


r 
r 

J" dx hi 


( k^q 

- e^- ld\ sin(kqx) sin(k^/x), (A3) 

2 m ' 


dxA(x)sm(kqX)sm(kq^x), (A4) 

(v) sin(^^x) sin(^^/ x), i = x, y, z, 

(A5) 


Appendix B: Spin Rotation Matrices 

Here we show how to perform the spin rotations for the 
two triplet components fo and /i (the singlet amplitude is of 
course invariant under spin rotations). The problem simplifies 
if all one wishes is to align the spin quantization axis with the 
local magnetization direction: this affords easier physical in¬ 
terpretation of the results. The central quantity that we use to 
perform the desired rotations is the spin transformation matrix 
T in particle-hole space. The quasiparticle amplitudes trans¬ 
form as. 


X(x) = mx). 


In our notation the matrix T can be written as: 


T = 


Ji 0 
0 S ’ 


(Bl) 


(B2) 


where 6^ is the kinetic energy on the y - z plane. Our nu¬ 
merical procedure for calculating the supercurrent involves 
initially assuming a constant amplitude form for the pair po¬ 
tential in each S layer, but with a total phase difference A(f 
(0, A(f at each S region). We then expand the pair potential 
via Eq. (A4). Similarly the exchange field and free particle 


where the submatrices and B are trigonometric functions 
solely of the angles that describe the local magnetization ori¬ 
entation. Expressing the orientation of the exchange fields in 
the regions Fi and F 2 in terms of the angles Of and 0/ intro¬ 
duced in Eq. (6) we can write Ji and !B as the following 2x2 
matrices: 
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cos (0//2) sin (O'l '^ + i sin (0//2) sin(^.) - cos(0//2) sin(^.) - i sin(0//2) sin(^^^) 
cos(0//2) sin(^p - i sin(0//2) sin(^^:^) cos(0//2) sin(^^^) - i sin(0//2) sin(^p 


(B3) 




cos ((pi/2) sin (O'l ') - i sin (0//2) sin(^.) cos(0//2) sin(^.) - i sin(0//2) sin(^^^) 
- cos(0//2) sin(^p - i sin(0//2) sin(^^:^) cos(0//2) sin(^^^) + / sin(0//2) sin(^p 


(B4) 


Here we have defined = Oijl ± 7 t/4. Using the spin rota¬ 
tion matrix 7“, we can transform the original BdG equations 
= 6n^n (Efi- (4)) by performing the unitary transforma¬ 
tion: 7^' = T1-IT~^ (of course we have T^T = 1). We then 
end up with the magnetization effectively along the new z axis 
and: 


f^o-h 0 0 A ^ 

0 + h A 0 

0 A* -^o-^h 0 

, A* 0 0 -^o-h, 


(B5) 


One of the benefits of working in this rotated coordinate 
system is that now the Hamiltonian matrix can be reduced to 
a smaller 2x2 size by using symmetry properties that now 
exist between the quasiparticle amplitudes and energies. 
As is the case under all unitary transformations, the eigen¬ 
values here are preserved, but the eigenvectors are modified 
in general according to Eq. (Bl). Thus for example, oper¬ 
ating on the wavefunctions using Eq. (Bl), and examining 
the terms involved in calculating the singlet pair correlations 
(Eq. (7)), it is easily shown that for a given set of quan¬ 
tum numbers n and position x, the following relation between 
the transformed (primed) and untransformed quantities holds: 

-h UniV^^. Thus, the terms that dic¬ 
tate the singlet pairing are invariant for any choice of quan¬ 


tization axis, transforming as scalars under spin rotations, as 
they should. 

The terms governing the triplet amplitudes on the other 
hand are generally not invariant under the spin-rotation. It 
is illuminating to see how both the equal-spin and different 
spin triplet correlations transform. The relevant particle-hole 
products in Eq. (9) that determine /o, upon the spin transfor¬ 
mations obey the following relationships: 

+ sin Oi [cos (pi - m„| v*^)+i sin (pi - m„| v* J]. 

Similarly the quasiparticle terms in the sum for /i (Eq. (10)) 
transform as: 

= si” J (B7) 

+ COS 0i [cos (pi - Un^vl^+i sin (pi [un[Vly - 

Thus the triplet amplitudes /o and f\ in the rotated system are 
linear combinations of the /i and /o in the original unprimed 
system (and vice versa). It is a simple matter to go from the 
rotated to the original system (and vice versa) by the route 
expressed in Eq. (Bl). 
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